Charge-2e Skyrmion condensate in a hidden order state 
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A higher angular momentum {£ — 2) d-density wave, a mixed triplet and a singlet, interestingly, 
admits skyrmionic textures. The Skyrmions carry charge 2e and can condense into a spin-singlet 
s-wave superconducting state. In addition, a charge current can be induced by a time-dependent 
inhomogeneous spin texture, leading to quantized charge pumping. The quantum phase transition 
between this mixed triplet d-density wave and skyrmionic superconducting condensate likely leads 
to deconfined quantum critical points. We suggest connections of this exotic state to electronic 
materials that are strongly correlated, such as the heavy fermion URu2Si2. At the very least, we 
provide a concrete example in which topological order and broken symmetry are intertwined, which 
can give rise to non-BCS superconductivity. 



I. INTRODUCTION 

It has become very much in vogue to argue that topo- 
logical aspects of condensed matter bear no relation to 
broken symmetries ji In a strict sense this need not be soi^ 
One can construct examples where a broken symmetry 
state has interesting topological properties and can even 
be protected by the broken symmetry itself. An inter- 
esting example of a mixed triplet d-density wave and its 
possible relevance to one of the many competing phases in 
the high temperature cuprate phase digram was recently 
demonstrated^ where it was found that the system ex- 
hibits quantized spin Hall effect even without any explicit 
spin-orbit coupling. In the present paper we show that 
the system exhibits charge 2e skyrmions, which can con- 
dense into a remarkable superconducting state. As we 
shall discuss, such a mixed triplet c?-dcnsity wave system 
and the resulting superconductivity is potentially rele- 
vant to the heavy fermion URu2Si2 with hidden order 

An early attempt at such a non-BCS mechanism of su- 
perconductivity was made by Wiegmanu)^ as an exten- 
sion of Frohlich mechanism to higher dimension. More 
recently, several interesting papers have led to discussions 
of superconductivity in single and bilayer graphenes. 
Grover and Scnthil^ have provided a mechanism in which 
electrons hopping on a honeycomb lattice can lead to 
a charge-2e skyrmionic condensate, possibly relevant to 
single layer graphene. To a certain degree we follow their 
formalism; see also the earlier work in Ref. of charge-e 
skyrmions in a quantum Hall ferromagnet. As to bilayer 
graphene, a charge-4e skyrmionic condensate has been 
suggested by Lu and Herbut^i and Mooni^ In the present 
paper we base our results, instead, on an unusual sponta- 
neously broken symmetry generated by electron-electron 
interaction, not by a given non-interacting band structure 
of a material, namely the mixed triplet-singlet density 
wave state of angular momentum i = 2, and we point out 
its possible implications to the mysterious hidden order 
state in URu2Si2, in particular to its superconductivity. 

It is appropriate to comment on what we mean by "hid- 
den order" . An order parameter can often be inferred 



from its macroscopic consequences in terms of certain 
generalized rigidities. Sometimes its direct microscopic 
signature is difficult to detect: a direct determination of 
superconducting order, a broken global U{1) gauge sym- 
metry, requires subtle Josephson effect, and even antifer- 
romagnetic order requires microscopic neutron scattering 
probes. Density wave states of higher angular momen- 
tum, such as the mixed triplet d-density wave, are even 
harder to detect. It does not lead to a net charge density 
wave or spin density wave to be detected by common s- 
wavc probes. It is further undetectable because it does 
not even break time-reversal invariance. A discussion of 
possible experimental detections of particle-hole conden- 
sates of higher angular momentum was given in Ref. [Tol . 
Thus, it is fair to conclude that the state we consider 
here is a good candidate for a hidden order. 

It is also necessary to remark on the realization of 
particle-hole condensates of higher angular momentum. 
An effective low energy theory of a strongly correlated 
system is bound to have a multitude of coupling con- 
stants, perhaps hierarchically arranged. In such cases, we 
can generally expect a phase diagram with a multitude 
of broken symmetry states. It is a profound mystery as 
to why non-trivial examples are so few and far between. 
A partial reason could be, as stated above, that these 
states are unresponsive to common s-wave probes em- 
ployed in condensed matter physics and therefore appear 
to be hiddeufii 

The next question is: are these low energy effective 
Hamiltonians contrived? If so, it would be of value to 
pursue them. However, simple Hartree-Fock analyses 
have shown that they certainly are not i^°i^^ an onsite 
repulsion U, a nearest neighbor interaction V, and an ex- 
change interaction J are sufficient in a single band model. 

The structure of this paper is as follows: in Sec. [Hi 
we construct the low energy effective action of the mixed 
triplet and singlet d-density wave system. In Sec. IIIIl 
we compute the charge and the spin of a skyrmion and 
verify that the skyrmions in this system are bosons, which 
can lead to a superconducting phase transition. In Sec. 
IIVI we compute the angular momentum of a skyrmion. 
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In Sec. |Vl wc study the charge pumping due to a time- 
dependent inhomogeneous spin texture that is interesting 
in its own right. In Sec. lVII we discuss mainly the problem 
of URu2Si2. In the appendices, the derivation of the non- 
linear (T-model and the details of computing the Cliern- 
Simons coefficients and charge pumping are provided. 



II. EFFECTIVE ACTION 

In the momentum space the mixed triplet and singlet 
d-density wave order parameter is (c and are fermonic 
annihilation and creation operators, respectively, Q — 
(7r,7r), and the lattice constant is set to unity) 



(1) 



where A'^ is a unit vector, a are the Pauli matrices acting 
on spin indices, and the form factors 



Wk 
Afc 



(cos k 

2 ^ 



cos ky ) , 
Aq sin kx sin ky, 



(2) 
(3) 



correspond to the d^i^yi and d^y density wave, respec- 
tivelyi^. It is not necessary that d^y and d^i_yi transi- 
tions be close to each other, nor are they required to be 
close in energy^. 

If we choose the spin quantization axis to be z, the up 
spins represent circulating spin currents corresponding to 
the order parameter d+id and the down spins to d—id (in 
an abbreviated notation). So, there are net circulating 
spin currents alternating from one plaquette to the next 
but no circulating charge currents. By the choice of the 
quantization axis we have explicitly broken SU{2), but 



not J7(l), and the coset space of the order parameter 
5^ = SU{2)/U{1). Such a state can admit Skyrmions 
in two dimensions (ignoring the possibility of hedgehog 
configurations in (2 -f 1) dimensions (cf. below). 
The Hamiltonian is 



E 



Ac 



_ (4) 

where the summation is over the reduced Brillouin Zone 
(RBZ) bounded by ky zL kx = ±7r, the spinor is ipl, ^ = 



and ek = —2t{coskri. + cos fey); addition 



of longer ranged hopping will not change our conclu- 
sions^. Here {i = x,y,z) are Pauli matrices acting 
on the two-component spinor. It is not necessary but 
convenient to construct a low energy effective field the- 
ory. For this we expand around the points Ki = (f , f ) 
and K2 = (— ^, ^), what would have been the two dis- 
tinct nodal points in the absence of the dxy term, and 
A'3 = (0, tt), what would have been the nodal point in 
the absence of the dx2-y2 term. This allows us to de- 
velop an effective low energy theory by separating the 
fast modes from the slow modes. After that we make a 
sequence of transformations for simplicity: (1) transform 
the Hamiltonian to the real space, which allows us to for- 
mulate the skyrmion problem; (2) perform a 7r/2 rotation 
along the r^-direction, which allows us to match to the 
notation of Ref. [l3 for the convenience of the reader; (3) 
label iJjKi+q,a by ipia, since Kt is now a redundant no- 
tation; (4) construct the imaginary time effective action, 
with the definition "0 = —iip^r^'. Finally, after suppress- 
ing the spin indices, and with the definitions 7*^ = r^. 



and 7^ 



-T^, we obtain the effective action in 



a more compact notation: 



3=1,2 



Wn 

-ij°dr - 2^^7"(^7,a, + dy) + i^{a ■ N)^y{-T^,dx + dy) + njjA„ 



+ / d^xips -11° dr -Wq{u- N)jy V: 



(5) 



where 771 = 1 and rj2 = —I. 



III. THE CHARGE AND SPIN OF A 
SKYRMION 

Wc will compute the charge of the skyrmions in the 
system by following Grover and Scnthil's adiabatic ar- 
gument 1^ First, consider the action around Ki = (f, f) 
when the order parameter is uniform (say, N = z). The 
results for K2 = (^f i f ) a-nd = (0,7r) follow identi- 



cally. In our previous paper we have shown that in this 
case the non-trivial topology leads to a quantized spin 
Hall conductance in iadx^-y^ + dajy-density wave stated 
as long as the system is fully gapped. The spin quantum 
Hall eflfect implies that the external gauge fields A'^ and 
A" couple to charge and spin currents, respectively. In 
the presence of these external gauge fields, we add mini- 
mal coupling in the action by 

^d^=p^^p^ + A';^ + ^A'^^. (6) 
Then the action is 
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i^^'dr + + ^K) - 2itr{d. + dy) + 2t7-(^^ + -Al +AI + —A^ 



.,_£^-^.(_0, + dy) -fa^^y{-Al -A^ + Al + -A^ + zAq 



(7) 



where we set e ~ h = 1. The non- vanishing transverse spin conductance imphes that the low energy efFective action 
for the gauge fields is given by 



and the charge current is induced by the spin gauge field 

1 



Ztt 



(8) 



(9) 



Consider now a static configuration of the N field with points in space such that U^{a ■ N)U = a^, and defining 
unit Pontryagin index in the polar coordinate {r,d): ip ^ Utp'-, and t/i = "ip'U^, we obtain 

N{r, 9) = [sin a{r) cos 9, sin a(r) sin 6, cos a{r)] (10) 

with the boundary conditions a(r = 0) = and a(r — >■ 
oo) = TT. Performing a unitary transformation at all 

I 



Wo 



i^°dr - 2itj^id., + dy) + i^a^^yi-d, + dy) + iAo 



To proceed, we write down the explicit form for U{r, 9), which is 



.Wo 



i-f°{U^drU) - 2it-i^{U^d^U + U^dyU) + i^cr^7«(-C/t9^i7 + U^dyU) 



(11) 



U{r,e) 



i(r) 

"('•L ie 



a(r) —iff 

-sm— j^e 

a(r) 



(12) 



In the far field hmit, Wd^U = (^^fs^)cr^ and WdyU = (^-^)cr^ substituting into Eq. (HH and introducing 

ff, = -iWd^U, we get 



d^x ^[ 



-ij^dr - 2zt7"(a, + dy) + i^a'-^y{-d, + dy) + ^Ac 



2<7"(/. + /,) + -^aV(/.-/y) 



(13) 



Equating the above equation and Eq. ([7]), we obtain in 
the far field limit 



A- -A^-Q- A'- 2sing 2COS0 



In other words, the process of tuning the order parameter 
from cr^ to (7 • N{r, 9) is equivalent to adding an external 
spin gauge field 



The total flux of this gauge field is clearly 47r. Suppose 
we adiabatically construct the Skyrmion configuration 
iV(r, 9) from the ground state z in a very long time period 
Tp — ^ cx). During the process, we effectively thread a spin 
gauge fiux of An. The transverse spin Hall conductance 
implies that a radial current will be induced by the 47r 
spin gauge flux of A'^{t), which is now time-dependent: 
A^(t = 0) = and l"(i = r^) = A", that is, 



- 2sin6l 2cos6' 2 

A'* = X H y -9. 

r r r 



(14) 



frit) 



27r 



dtAm- 



(15) 
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As a result, charge will be transferred from the center to 
the boundary, and the total charge transferred is 



IV. THE ANGULAR MOMENTUM OF A 
SKYRMION 



dt 



2tt 



-2. 



(16) 



Therefore, after restoring the unit of charge to e, we ob- 
tain a Skyrmion with charge 2e; its spin is 0. 

It is important to verify the adiabatic result by a dif- 
ferent method. This can be done by a computation of the 
Chern number The charge and spin of the skyrmions 
are associated with the coefficients of the Chern-Simons 
terms by the following relations: Qskyrmion = €26 and 
'S'skyrmion = C*! | , where Ci and C2 are 



Ci 



247r2 



Tr 



dk.j 



(17) 



d*k{a-z)G- 



dka 



dki, dk\ 



(18) 



where G is the matrix Geen's function and the trace is 
taken over the spin index a and other discrete indices. 

If the Green's function matrix is diagonal in the spin in- 
dex, then the Chern-Simons coefficients for up and down 
spins can be computed separately. 



Tr 



dkx 

(19) 

and Ci = A/'(G|) + A/'(G^), C2 = Af{G^) - A/'(G^). Fur- 
thermore, it can be shown (see Appendix IB|) that for 



g: 



iuji 



(20) 



with hg- being the Anderson's pseudospin vector— of the 
Hamiltonian, the Chern-Simons coefficient for spin a can 
be written as 



where = /\h„\ is the unit vector of h„. Here Ci and 
C2 are the total Chern number and the spin Chern num- 
ber A/'spin defined in our previous paper, respectively.— For 
i(Tdr^2_y2 + dxy system, we have /lo- = (Afc, — crVFfe, e^). 
Exphcitly, Gi = -1 -M = and G2 = -1 - 1 = -2; thus 
the results are the same as above. 

Because a Skyrmion in the system carries integer 
spin, it obeys bosonic statistics and may undergo Bose- 
Einstein condensate. As a result, the charge- 2 e Skyrmion 
condensate will lead to a superconducting phase transi- 
tion. But what about its orbital angular momentum? 
In the following section, we will prove that it is zero re- 
sulting in a s-wave singlet state. This is a bit surprising 
given the original d-wave form factor. 



To compute the angular momentum carried by a 
skyrmion in the system, we consider the angular momen- 
tum density due to the electromagnetic field. For a static 
spin texture it is clearly zero, because E ~ 0. For a time 
dependent texture it is little harder to prove. Consider, 

Nx{r,9,t) = ama{r,t) cos I3{e,t), 
Ny{r,0,t) = sin a{r,t) sin /3{0,t), 
Nz{r, t) = cosa(r, i), 

where a(r, t) and /3(0, t) are smooth functions, and a(r, t) 
satisfies the boundary conditions a{r = 0,<) = and 
a{r 00, t) = TT, for any t, and ^^^l^^cx) = 

^"at *^ |r- »oo = in the far field limit. The unitary matrix 
is now time dependent. After a little algebra, we obtain 
the time-dependent gauge fields in the far field limit to 
be 



So, ^9,t) = At{9,t) 



-2 sin 6* (9/3(6*, t) 
r d9 ' 
2 cos 9 d/3{9,t) 



89 



(22) 
(23) 



~ ^ at 



and A%r,9,t) 



AJ(r, 9, t)x + At{r, 9, t)y = A^(r, 9, t)9, where 



A^eir,9,t) 



2 dpi9,t) 
r 89 ' 



(24) 



Therefore, the electric field will have a non-zero 
0— component, E = Eo9^ and the magnetic field will have 
a non-zero z— component, B = BzZ, where 

1 8Af{9,t) _ 8A'g{r,t) _ 4 8^ (3(9, t ) 
^ 898t 



Eh 



r 89 8t 
8AI{r,t) 2 8/3{9,t) 



(25) 
(26) 



8r 89 
As a result, the angular momentum density still vanishes, 

1 



Id 



47rc 



fx {Eg9 X Bzz) = 0. 



(27) 



It is possible that superconductivity with non-zero an- 
gular momentum may be realized when the interaction 
between skyrmions is included, but we do not know how 
to prove it. It would be interesting to explore what other 
kinds of quantum numbers are carried by the topological 
textures in the model we have studied. 



V. QUANTIZED CHARGE PUMPING 

In Sec llin we considered a static spin texture and ob- 
tained charge-2e skyrmions in the system. If we consider 
a time-dependent spin texture, which has a slow varia- 
tion in one spatial direction, say, y, and is uniform in the 
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other. X, charge will be pumped from one side of the sys- 
tem to the other along xJ^ This charge pumping effect 
can be understood from the effective gauge action, which 
is 



(28) 



where the integral is over the real time, t, instead of the 
imaginary time, r. Therefore, the charge current induced 
by the spin gauge field will be 



J, 



C2 
47r' 

Stt' 



(29) 



where we define the spin gauge flux F^^^^ = d^Af, — d^A^j^. 
After some straightforward algebra (see Appendix[C]) , the 
spin gauge flux can be written in terms of the iV- vector. 



(30) 



As a result, even in absence of the external electromag- 
netic field, a charge current may be induced by a time- 
dependent inhomogeneous spin texture because 

C2 



Stt 



E^'^^N -[{d.N) X {dxN)]. 



(31) 



To demonstrate the charge response induced by the 
spin texture, we consider the following configuration with 
unit Pontryagin index, 

N{y, t) = [sin6'(i) cos (/)(?/), sin 0(<) su\(f){y), cos9{t)] , 

(32) 

where 9{t) and 4'{y) are smooth functions of t and y, 
respectively, with boundary conditions 9{t = 0) = 0, 
9[t = Tp) = TT, and (f){y — ±c») = ±7r. Therefore, we 
have an induced charge current along the i-direction, 



Jx 



"""^N ■ [{d^N) X (dxN)] 



Stt 

■ [{dyN) X (dtN)] 



(33) 



Interestingly, we can show that the pumped charge is 
quantized. 



Qpumpod = / dt dy 
Jo J -00 



^ I dt I dyN- [{dyN) X (a* TV)] 
C2 



^ , d9 d(j)N- [{dgN) X {d^N)] 

= C2, (34) 

where we have used that, for the spin texture with unit 
Pontryagin index, 

/ d9 I d(j)N- [{deN) X {d^N)] = 47r. (35) 

Jo J-TT 



After restoring the unit of charge, we have Qpumpcd = 
C2 e. So far we have considered the spin texture with 
unit Pontryagin index. If the spin texture is generalized 
to a general Pontryagin index, Np, then the pumped 
charge will be Qpumpod = C'2Np e. 

How could we observe this charge pumping experimen- 
tally? We need to control the direction of the iV-vcctor 
so that it can be the time-dependent inhomogeneous spin 
texture discussed above. In topological chiral magnets,— 
the TV- vector is the net ferromagnetic moment, which 
aligns along the external magnetic field, so one can ap- 
ply a time-dependent magnetic field H{t) = H(t)x cou- 
pling to the A'^-vector and control the magnitude of x- 
component of N. 

In the mixed triplet d-density wave, however, the sit- 
uation is more complicated. In the presence of an ex- 
ternal magnetic field, there will be a spin flop transition 
and the iV-vector will lie in the plane perpendicular to 
the external field.— In other words, we cannot fully con- 
trol the direction of N with a time-dependent magnetic 
field. Therefore, it would be a challenge to measure the 
pumped charges in the system. 

Nevertheless, the charge pumping effect provides, at 
least, a different conceptual approach to probe the topo- 
logical properties of the system in addition to the quan- 
tized spin Hall conductance. For the quantum spin Hall 
effect, a spin current is induced by the external electric 
field,?', whereas for the charge pumping effect, a charge 
current will be induced by the spin texture. It would, of 
course, be interesting if one can manipulate the A^-vector 
experimentally because the charge current is easier to de- 
tect than the spin current. 



VI. DISCUSSION AND APPLICATION TO THE 
HIDDEN ORDER STATE IN URu2Si2 

There are two points that we have glossed over. The 
first is rather simple: in the ordered phase at T = 0, there 
arc also Goldstonc modes that can be easily seen by inte- 
grating out the fermions resulting in a non-linear cr-model 
involving A^, the form of which is entirely determined by 
symmetry. These do not lead to any interesting physics, 
such as charge-2e skyrmions that condense into a super- 
conducting state. At finite temperatures they could lead 
to a renormalizcd classical behavior.— The second point 
is more subtle: we have assumed that the hedgehog con- 
figurations are absent. This would require, as pointed out 
by Grover and Senthil^^, that the energy of the Skyrmion 
(especially in the limit Aq — >■ 0) is smaller than indi- 
vidual pairs of electrons, a question that is likely to be 
model dependent. If this assumption is correct, however, 
the transition from the the mixed d-dcnsity wave state 
to the superconducting state will correspond to a decon- 
fined quantum critical point, which otherwise would have 
been a first order transition, as in Landau theoryJi 

We suggest that the superconducting phase driven by 
the skyrmion condensate may be realized in the URu2Si2, 
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which hosts an exotic hidden order (HO) phase, with bro- 
ken translatenal symmetry below T^q ~ VI .^K and a 
superconducting phase below Tc ~ Recently, Fu- 

jimotoi^ proposed a triplet d-dcnsity wave with the or- 



der parameter (c[ i „Cfc+Qo, 2,/3) 
i{Ai sin 



d(k) ■ Sap with d{k) = 
— sin fc^, 0,0) to describe this state^^ here 

1 and 2 refer to two different bands and Qo = (0, 0^) is 
the nesting vector; even the earlier work in Ref. |l9j in- 
volving circulating spin current is not entirely unrelated. 
The order parameter considered in Ref. [l^ is different but 
a close cousin of the order parameter considered in our 
work; the circulating staggered spin currents in Ref. [l^ 
lie on the diagonal planes instead and the crucial dxy part 
is missing there. That the currents are in the diagonal 
planes instead of being square planar is conceptually not 
important, but is necessary to explain the nematicity ob- 



served in the experiments!^ We now discuss the role of 
spin-orbit coupling before making our final comments. 



Spin-orbit coupling 

It will be shown below that the order of magni- 
tude of the spin-orbit energy i^so ~ [(-^ • 5)^ — 
l](A2/Wl(H^o/W')2[l-hCl(W^o/W^)^], correcting a mistake 
in Ref. [ijj. Here Ao is the strength of the spin-orbit cou- 
pling, given by 



(36) 



where A (A:) = (Ao/\/2)[a;sinfcy— ysinfc^:]. In the presence 
of spin-orbit coupling, the Hamiltonian is 



■Htotal = "H + 'Kso 

eu Afc + iN.Wu 

2^ "^fc I A,(fc) + iKy{k) lWk{N, + iNy) 

-iWk{Na, + iNy) ~A4k)~iAy{k) 



K{k)~iKy{k) iWk{Nx-iNy) 

-iWk{Nx - iNy) -A^(fc) + iAy{k) 



iN.Wk 



(37) 



where ^'j^. is the four-component spinor (c]^, ^, c^f.+Q f ''fc J.' '^k+Q 4.)' absence of spin-orbit coupling, the eigenvalues 

are ±Eok with Eok = V^fc + ^k + ^1- other hand when spin-orbit coupling is present, the eigenvalues of the 

upper and lower bands now become Aup,± = Ek.±, Aiow.i = ~Ek.±, respectively, where 



E, 



k.± 



A? 



{4 + Wl)Kl-W^,{N-KkY 



(38) 



with A^ 



lA, 



A^(fc) 4- A^(fc). When the dxy component is absent, = 0, and the results of Ref. [l^ are 



recovered. Consider the following two cases separately. 



iV 



Since the chemical potential is at the mid-gap, we can 
focus on the lower bands. When N = z,we have A^- A^ = 
and 



Al±2[Ei,Alr 
= -EokT\Ak\ 



(39) 



I 

ergy will be 



n- 



E, 



Ok 



\Au 



) - 2{-Eok)] 

Eok + |Afe|) + 2Eok 



(40) 



2. N ±z 



When N lies in xy-plane, we have N ■ A^ = | A^ | cos (j)k , 
where (j)k is the angle between N and Afc, and 



Assuming that Aq ^ Wo , Aq <C W with the electronic 
bandwidth W = 8t, the change in the ground state en- 



cos (pk 



iV- Afc ^ iV,A,(fc)+A^^A^(fc) 
lAfcl y^A2(fc)+A2(fc) 



(41) 
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The eigenvalues of the lower bands are now 



A|±2[£;„\A2-M/2A|cos2 



-^„,^(l-__J)|A.|--^^(l + 0(^)), 



^Ofe 



where we have used cos^ ( 



0(1). Notice that the signs of the second order terms for Aj^^ ^ and Aj^^^ _ 



(42) 



are both 



negative, leading to the net change in the ground state energy, which is opposite to the N ~ z case. Assuming that 
Aq <C Wo, Aq <C W, the change in the ground state energy per lattice site will be 



so 



y 

low, 



+ ^low.- 



2{-E, 



Ok) 



A2 

'W 



^Ofc 

Wo 
W 



l + 0(-f) 

^Ofe . 



l + O 



w 



< 0, 



(43) 



Therefore, iV-vector should lie in the xy-plane in the presence of spin-orbit interaction and the result stated above 
follows. 

I 



As large as the spin-orbit coupling may be for U atoms, 
i?so is still a small energy scale. However, if other 
anisotropies are absent, the order parameter would be in 
the AK-plane, resulting in vortices; exchange anisotropy 
can also result in an easy-axis anisotropy, in which case 
spin textures could be Ising domain walls that can trap 
electrons. Although skyrmions are finite energy solu- 
tions, vortices cost infinite energy unless they are bound 
in pairs. We speculate that charge 2e-skyrmionic con- 
densation is a more likely scenario, but the crossover in 
the texture is an interesting topic for further research. 
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The following remarks about URu2Si2 are relevant: in 
both magnetic field-temperature {H — T) and pressure- 
temperature (P — T) phase diagrams, the superconduct- 
ing phase is enclosed within the HO phasei^ It implies 
that the superconducting phase is closely related to the 
HO phase, and is probably induced by it. Throughout 
our calculation, ignoring of course skyrmions, we have 
assumed that the system is half-filled. The lower band is 
filled and the upper band is empty, and the topological 
invariant is quantized. If this is not the case, then there 
will be no quantized spin Hall conductance, but an in- 
duced superconducting phase from charge 2e-skyrmionic 
condensation; doping will result in conducting mid-gap 
states, as in polyacetylenei^ Of course, such a topolog- 
ical superconducting phase is very sensitive to disorder. 
Indeed, this may be supported by the destruction of the 
HO and SC phases with 4% Rh substitution on the Ru 
site4. To summarize, we can find a rationale for a hidden 
order phase enclosing a superconducting phase at lower 
temperatures. 



Appendix A: Derivation of the non-linear (T-model 

To derive the non-linear cr-model, we compute the ef- 
fective action by integrating out fermions. We start with 
the action S ~ where 

S, ^ J d'x^, [G-i] Vj, (Al) 

with GJ^ = Goj 

For j ~ 1,2, we have 

Goj = -ia'^T'dr - 2ita°Ty{iyd^ + dy), (A2) 
EE 1,7, AoaV" - i^{a ■ N)T%-7j,d, + dy)iA3) 
and for j = 3, we have 

Ggj = -lCT°T'dr, (A4) 

E3EE Wo(CT-7V)r^. (A5) 
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The effective action will be Seff = '^j.eff with where the fermion operators can be integrated out easily 

since the Hamiltonian has only bilinear fermion opera- 
tor terms. Using the mathematical identity Indet = 
trln^ with tr being the trace, we have 



S. 



In 



DtpjDtpje 



-ln[det|G-i| 



(A6) 



%e//=-trlnGoj[l + Go,,E,] 

= -trlnGoj - tr [Go,,S,] + ^tr [Go,,S,Go,,S,] 



(A7) 



r 



where we have used ln(l + x) = x — ^ + ■ ■ ■ . particles and the first order term vanishes, so our goal is 

The zeroth order term is the effective action for free to compute the second order terms: 

I 



dr / dr' d x d x'Tr [Gqj (x, r; x', r')Ej(a;', r')Goj (x', r'; x, r)Ej (x, r)] 



i ^ Tr [Go,,(fc)S,(g)Go,,(fc + q)^,{~q) 



(A8) 



k,q 



r 



where k = {ko,kx,ky), q = (qa, qx, qy), and Ga_j{k) can Putting all together, taking long wavelength limit (q — > 
be obtained by inverting Ea. (|A2[) and Eq. (jA4[) . 0) and keeping only terms up to the second order deriva- 

tive, we have, for j = 1,2, 



c(2) ^ 2 V - 

J,eff- ^ kl+Afii'q.kx + ky) 

k,q 



-Al + {^ri-,yqx + qynNg-N.g) 



(A9) 



where terms which are odd in k and q are dropped. 

Using the relation J2qfqf-q = J dTd'^x\f{x,T)\'^ , we 
obtain 



c(2) , q{2) 



1 

.91 



drd X 



dxN 



dyN 



(AlO) 

where the constant terms are dropped, (AT, Y) is the co- 
ordinate after a 7r/4 rotation, and 



-Wo' 



51 ^2ikl + At^ikx + ky)^) 
Similarly, for j ~ 3, we obtain 



(All) 



k.q 

drd^x 

93 



drN 



(A12) 



where 



(A13) 



Therefore, we obtain the non-linear sigma model. 



5*6// ~ — / drd'x 



(A14) 



where the constant terms and higher order terms are 
dropped, and it is rescaled in order to obtain a famil- 
iar form. 
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Appendix B: Chern-Simons coefficients 



We start by taking (/i, ly, A) to be (0, x, y), and obtain 



In this appendix we are going to prove that 



247r2 



Tr 



dk,, 



4n dkx dky 



G, 



dG- 



(Bl) and 



1 



1 



[iLoI + f • ha) ■ (il) 

{^ui + iT-K), (B2) 



G, 



dG- 



1 



dk^ [i^Y _ \h^\2 
-1 



iuji + f ■ ha- 



^-^■dk-J 



it dha f , , dha , - dha 



(B3) 



where we have used the matrix identity (f ■ a){T ■ b) = (a ■ h)I + ir ■ (a x b). Similarly. 

-1 



Ga- 



dky (iu^y - 



dha.f . dh 



^ ^ dha_^ 

dky dky 



(B4) 



Therefore, 



dha., dh. 



dha , , dh„ 



dh„ 



dkx dk 



^ ■ i^'ai 'r ha- X ——) 



dky 



dky 



(B5) 



Since we are going to multiply it with the antisym- (x o y) will vanish. Therefore, only the last term in the 
metric tensor e^^A, the terms which are symmetric under bracket contributes, 



- I ^'^cr , r dha. 



dkx 
dha 



dkx ' 



^ I dha , - dha. 



dk 



dha 



d kj^ d ky 



- iohai 



dha 



y 
dha 



dkf dkqi 



dky 



dha .J- dha. 

■ i^^T7—(ha ■ ) 



dky 



dkx 



J- ,dha dha. , dha ,r: dha. , dha dha. 7' 



dk^ dky 



dkx 



dky 



d krf^ dky 



(B6) 



where we used the following mathematical identities: Therefore, after combining with e^xy and taking the 

^ ^ ^ trace, we have 

ax {b X c) ^ b{a ■ c) - c{a ■ b) , 

{a X b) X {a X c) ^ {a ■ (b X c))a. 
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■ ac-i 

Lrrr — 1^ ^tr—^, <-J( 



9w 



dk,, 



(M2-|/l,P)3 



-Tr < —iujT ■ 



-{r-K) 



,dha dh, 



dhcr dh„ J' 

■ dk: dk-J''" 



dh„ dh„ 



(M2-|/^,|2)3 
-2 



dha- . 
d kx ^ ky 



— (h ^^"^ dhcr 

((iw)2-|/v|2)2 ""dk^"^^' 



(B7) 



where we have used the fact that Pauh matrices are trace- 
less, so the only contribution will be the term propor- 
tional to /. 

We have six non-zero terms because of the e^uX tensor, 

so 



N{G„) ^ - 



2 • 6 



247r2 



d^k 



1 



{{iujy - \K\^Y dk^ dk 



(fk 1 ,^ dhc dhcr-. 



4^ l/i.P 



fk- 

47r 



dha- dhcr 
dhjx dkfj 



(B8) 



where the energy integral was done by computing the 
residue of the second order pole. 



where 



'iWdf^U. 



Therefore, we can write the spin gauge field in terms 
of the unitary matrix. 



and we have 



^^ = Tr 



2 



-iTr [a'U^'df.U] 



(C2) 



Appendix C: Spin gauge flux in terms of N 

In the main text, we obtain the spin gauge field to be 



fu 



2 



(CI) 



J 



= -zTr [a'{d^U^)id,U) - a' id,U^)id^U)] (C3) 



Assume that the spin texture has the general form: 

N{x,t) = [sin9{x,t) cos(j){x,t),sm9{x,t) sm(t){x,t),cos9{x,t)] , 
where 0(x, t) and t) can be any function of position and time. Then, we have the unitary matrix 



U{x,t) 



cos 



sin^e^-^^-^'*) 



cos 



e{x,t) 



(C4) 



(C5) 



and 



d^uHx-.t) 



d,U{x,t) 



\ sin 1 9^ 6 



e'^-icos|a^ 



2 v.'-'a 2^^^ 2 sin 2^^^ 



e-**(icosfa^6l-isinf9,,^ 
— ^ sin Tfdp^O 



— I sin ^duO 
e*'^(icosfa^6'-|-isinfa^(; 



e-*'^(-i cos \dj + i sin f 
— ^ sin ^dyO 



(C6) 



(C7) 
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where we have suppressed the arguments of 9{x,t) and 4>{x^t). 

Therefore, we can calculate the product of the last two matrices, and express the spin gauge flux as 



Ff. 



s'mO {d^9di,<i> — d,y9df_, 



X 2 



In addition, we can also write N ■ {dfj,N x d,yN) in terms of 6{x, t) and (j){x, t), 



N ■ (d^N X d^N) 

sin 6{x, t) cos t) 

[cos 9{x, t) cos (f>{x, t)df^9{x, t) 

— sm9{x, t) sin (j){x, t)d^(j){x^ t)] 

[cos 9{x, t) cos (j)[x, t)d^9{x, t) 

— sin6'(x, t) sin t)d^4)(x^ t)] 



sin t) sin t) 

[cos 9{x, t) sin t)dfj_9{x, t) 
'r sm9{x, t) cos (j){x, t)df^(j){x, t)] 

[cos 9{x, t) sin (/>(x, t)d^9{x, t) 
\- sm9{x, t) cos t)dy(f){x, t)\ 



cos9{x, t) 
ain9{x, t)d^9{x, t) 

suy9{x, t)dy9{x, t) 



= s\tl9 {d^9du4' - d^9dfj_(j)) , 
where, again, we suppressed the arguments of 9{x,t) and (j}{x^t). Finally, we obtain 



(C8) 



(C9) 
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